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ABSTRACT 


In this study, seismic analysis of soil-foundation interaction under a bridge pier are studied with different earthquake 
excitations. In 2016, August 25, a magnitude of 6.8 happened near Bagan region. Pakokku Bridge, the longest of the bridges 
over the Irrawaddy, is situated on 37.8 miles from the epicentre of 2016 Chauk earthquake. That is why the safety performance 
of long-span Bridge (Pakokku Bridge) especially for the safety of the foundation system subjected to soil-foundation interaction 
is necessary to investigate for unexpected future seismic excitation. Firstly super structural loadings on the pile cap are 
estimated by using STAAD PRO V8i. And then, p-y curves are determined by Reese (1974) method for the static and National 
Cooperative Highway Research Program (NCHRP) for dynamic conditions. Based on the development of p-y curves, theoretical 
ultimate soil resistance per and ped due to wedge and flow failure are determined to produce critical depth xcr. After that, finite 
element software ABAQUS is used for the analysis of soil-foundation interaction under a bridge pier in static condition. And 
then, the behaviour of soil foundation interaction under a bridge pier is carried out due to Chauk earthquake. In this study, the 
behaviour of soil-foundation interaction such as deflections and settlements are produced. According to the analysis results in 
static condition, it is found that the vertical and horizontal displacements at the pile tip are 2.28mm and 0.14mm respectively. 
In dynamic condition, the vertical and horizontal displacements at the pile tip are 3mm and 2.94mm are found at 0.12g. After 
that, maximum ground acceleration of 0.5g is 14.5mm and 2.94mm in horizontal and vertical displacement of pile tip. 
Maximum shear stress and strain are found out the base of the pile cap. Finally it is found that the soil-foundation interaction 
under a bridge pier during earthquake motions presented in this study is reliable and reasonable with the limitation of 
AASHTO Standard Specifications for Highway Bridges. 
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I. INTRODUCTION 

Pakokku Bridge is a rail and road bridge across the 
Irrawaddy River in Myanmar. The bridge is part of the India- 
Myanmar-Thailand Trilateral Highway and is the longest 
bridge in Myanmar. This bridge is connected to Pakokku city 
and the administrative district of Nyaung-U and Mindat. The 
construction was begun at 2009 and it was completed in 
2012. Pakokku Bridge is 3.4 km long and 14.8 m in width. In 
2016, a magnitude of 6.8 Chauk earthquake happened near 
Bagan region. Pakokku Bridge is situated on 37.8 miles from 
the epicentre of 2016 Chauk earthquake. Therefore, the 
safety performance of Pakokku Bridge especially for the 
safety of the foundation system subjected to soil-foundation 
interaction is necessary to evaluate due to seismic excitation. 


The research described in this paper presents a numerical 
investigation of composite soil-foundation system. Firstly, 
checking the bearing capacity of Case Study Bridge is based 
on Shamsher Prakash [2] and compare pier loading and pile 
group capacity. After that, p-y curves are determined using 
the subgrade reaction approach according to Reese and 
Matlock [6]. The soil-foundation system is modelled and 
analyzed using Abaqus/Cae [4]. 

II. RESEARCH METHODOLOGY 

The research study included structural investigation, soil 
investigation, modelling and analysis for the selected case 
study. The flow chart of this study is shown in Fig. 1. 


When lateral loads are applied on a pile, lateral deflection of 
the pile depends on the soil resistance, in turn, depends on 
the pile deflection and this dependence is known as soil-pile 
interaction. Pile foundation is one of the most common and 
important sub-systems of the bridges. Such a foundation is 
commonly chosen when the construction site is a weak or 
saturated sandy ground. Bridge-pier foundation is important 
to simultaneously consider the interaction behaviour of a 
soil-group pile-bridge pier system when designing such a 
foundation. Bridge substructure is a very important part of a 
bridge as it safely transfers the loads from the 
superstructure to the earth [1]. The subgrade reaction 
approach provides the simplest solution for the pile-soil 
interaction problem. The subgrade reaction has been widely 
accepted in the analysis of soil-structure interaction 
problems. The p-y approach is another method for handling 
pile-soil interaction. Moreover, the finite element method 
(FEM) is the most powerful tool in modeling soil-structure 
interaction. 



Fig. 1 Flow chart of the methodology 
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III. STRUCTURAL INVESTIGATION 

The selected case study, Pakokku Bridge is the longest 
bridge over the Irrawaddy River in Myanmar. It is composed 
of 6 main spans and 20 numbers of piers under the bridge. 
The superstructure of the bridge is composed of steel truss 
members, concrete deck slab, and I-beam girders. The 
behaviour of soil foundation interaction are analysed for the 
pier number 5 under the bridge of the middle span. The 
photo of Pakokku Bridge and case study of pier (PR 5) 
location, sectional elevation and plan view are shown in Fig. 
2 . 
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Fig. 2 (a) Pakokku bridge, (b) Study portion of Pakokku 
Bridge, (c) Sectional elevation of longitudinal and transverse 
of case study and (d) Plan view of case study (PR5) 

As shown in Figs. 2 (c) and (d), pier under the bridge are 
composed of two column bents and supporting foundation 
with eighteen piles group.. So there are two loading points, 
pier 1 and pier 2, on the pile cap. Table I and II present the 
information of pile and pile cap details, and material 
properties. 

Superstructural loadings on the pile cap are determined 
according to AASHTO. The factored loading from the 
superstructure on the pile cap are presented in Table III. 


TABLE I Pile and Pile Cap Details 


Pier No. 

No. of Pile 

Pile Cap Dimension 
(m x m x m) 

Diameter of Pile (m) 

Pile Length 
(m) 

5 

18 

24.6x11.8x2.75 

1.67 

18 


TABLE II Material Properties 



Pier 

Pile 

Design strength of concrete (MPa) 

24 

29 

Yield strength of re-bar (MPa) 

290 

290 


TABLE III Loading Of Structure On The Pile Cap 


Pier 

F x 

(kN) 

Fy 

(kN) 

F z 

(kN) 

M x 

(kN-m) 

My 

(kN-m) 

M z 

(kN-m) 

1 

0 

2.4xl0 4 

8.9 

18 

9.3 

1127 

2 

0 

2.7xl0 4 

0 

85.2 

8.3 

1537 


I. SOIL INVESTIGATION 

Standard penetration test SPT N values, friction angle ((])), the compressive strength of soil (c) and unit weight of soil (y) are 
important parameters to analyze the behaviour of soil-pile interaction. The properties of soil under the pier 5 are shown in 
Table IV. As shown in Table IV, the cohesion of soil V, for every layer is zero and the internal friction angles of soil are ranging 
from 19.5° to 30.5°. 
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TABLE IV Properties of Soil 


Depth z, (m) 

N 

Y(kg/m 3 ) 

d)fdegree) 

c(KN/m 2 ) 

3 

10 

1853.4 

19.5 

0 

6 

19 

1866.2 

19.5 

0 

9 

30 

1912.6 

19.5 

0 

12 

23 

1883.8 

19.5 

0 

15 

25 

1845.3 

21.3 

0 

18 

22 

1853.4 

21.3 

0 

21 

25 

1895.0 

21.3 

0 

24 

30 

1911.0 

29.2 

0 

27 

29 

1907.8 

29.2 

0 

30 

29 

1903.0 

29.2 

0 

33 

55 

1946.3 

29.2 

0 

36 

75 

1959.1 

30.4 

0 

39 

88 

1965.5 

30.4 

0 

42 

118 

1984.7 

30.4 

0 

45 

124 

2008.7 

30.4 

0 


II. CALCULATION OF THE BEARING CAPACITY OF 
BORED PILE 

The bearing capacity of piles is, therefore, estimated based 
on initial strength and deformation characteristics of the soil. 
The effect of changed soil conditions are reflected in the 
nondimensional empirical coefficient N q and mobilized shaft 
friction, f s , in cohesionless soils [2]. The final expression for 
ultimate load capacity, (Q v ) u it of a pile then becomes 
(Qv)uit = Qp +Qf (1) 

The ultimate end resistance Q p in tons of drilled or bored 
piles can be estimated by the following relationships 


For Sand, Q =-(o.4N/B)DfA p ^ — 


For cohesionless or nonplastic silt, 


NA P 


4(o 


1 L 


4N/BjD f A p < NA P 


( 2 ) 


( 3 ) 


The friction capacity of a pile 

The friction capacity of a pile can be estimated by using the 
following relationship: 

Q f = -^-(f s )(perimeter) (Embedment length) (4) 


evaluated form the basic procedure for developing p-y 
curves in the static condition (Reese et al., 1974) [6]. 

A. P-y Curves of Laterally Loaded Piles in Cohesionless 
Soils for Static Condition 

The p-y method models a laterally loaded pile as a soil- 
structure interaction problem because the lateral load 
applied results in lateral deflection of the pile, which causes 
reactions in the soil. A static equilibrium between the pile 
and soil must be obtained. The numerical solution of the 
problem requires a relationship between the pile deflections 
and the soil reactions. Static p-y curves with different depth 
of soil for transverse and longitudinal directions of the pile 
group are developed according to Reese and Matlock [6]. 

Soil resistance and deflection are determined by using the 
following equations. 


= yx 


K xtan^inp + _^ (B + Xtanptang) ~ 


_tan(p-^)cosa tan(p — </>) 

+ Yx[ko xtanp(tan<j)sin|3 - tana) - K,\ Ei] 

Pcd = K a Byx(tan 8 (3 -1)+ Ko Byx tan 4> tan 4 P 


(7) 

( 8 ) 


f s = — < ltsf (5) 

50 

where f s is the ultimate unit shaft friction in tons per square 
feet. 

(Q p )g = n Q P > (6) 

where (Q p )g is ultimate point load of a pile group, (Q p ) is the 
ultimate point load of a single pile and n is the number of 
piles [2]. 

The allowable bearing capacities of single pile and group 
piles of case study under the pier (PR5) are 3825 kN and 
68858 kN. The pier loading from superstructure, 51000 kN, 
shown in Table III is less than the allowable bearing capacity 
of group piles. Therefore, design of foundation system under 
(PR5) is satisfied. 

III. DEVELOPMENT OF P-Y CURVES FORSTATIC AND 
DYNAMIC CONDITION 

For the solution of the problem of a laterally loaded pile, it is 
necessary to predict a set of p-y curves. The p-y curves were 


Where ,a = ^§, p = 45 + a,and KA = tan 2 f 4 5-^-c|) 


Pm = BlP c 

(9) 

B 

ym_ 60 

[10] 

c£ 

< 

ll 

cu 

(11] 

_3B 
y " _ 80 

(12] 

m= Pu_Pm 

y u -y m 

(13] 

n = P ” 

m y m 

(14] 

p _ Pm 

~(yi 

(15] 

f c r /(n_1) 

yk= tJ 

(16] 

p = C y/n 

(17] 


@ IJTSRD | Unique Reference Paper ID - IJTSRD21339 | Volume - 3 | Issue - 2 | Jan-Feb 2019 


Page: 387 


































International Journal of Trend in Scientific Research and Development (IJTSRD) @ www.ijtsrd.com elSSN: 2456-6470 
where, 


Per, Pcd = theoretical ultimate soil resistance due to wedge 
and flow failure 

p c ^govern theoretical ultimate soil resistance 

p u = ultimate soil resistance 

y u = deflection at ultimate soil resistance 

p m = soil pressure at D/60 

y m = deflection at soil pressure at D/60 

p = establish initial straight line portion 

yk = deflection at establish initial straight line portion 

x = depth below the pile head 

y = unit weight of soil 

§ = angle of internal friction 

K 0 = coefficient of lateral earth pressure 

B = diameter of pile 

a, p, K a are necessary soil parameters for obtaining the 
value of x cr , at the interaction of p cr and p C d- The coefficients 
Bi and Ai are used for determining of p m and p u . The critical 
depth value x cr is obtained by plotting p cr and p C d with depth 
(x) on a common scale. 

Valuts of p cr and ft4 



Per 

Pcd 


Per and (k_\7m) 

Fig. 3 Critical depth 

Above the critical depth, p cr is taken for soil resistance and 
Pcd is taken below the critical depth for p-y curves [2]. The 
evaluation of critical depth to develop p-y curves are shown 
in Fig. 3. The developments of p-y curves for single pile of 
selected case study pile group are presented in Fig. 4. 
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Fig. 4 p-y curves for single pile 

Evaluations of p-y curves for pile-group are determined by 
applying P-multipliers P m [5]. The development of p-y curves 
for transverse direction and longitudinal direction of group 
piles are shown in Fig. 5. 


p-y curves at different depths fur pile gmup 
(Lcin^itudhial) 
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Fig. 5 p-y curves for (a) longitudinal direction and (b) 
transverse direction 

B. P-y Curves of Laterally Loaded Piles in Cohesionless 
Soils for Dynamic Condition 

The dynamic, single pile p-y curves were fit with an 
analytical expression, which appears to be valid for soft to 
stiff clay and loose to dense sand [5]. That expression is 
given below; 

Pd = * [« + + tea-, ] 


where, 


Pd 

Ps 

So 


y 

D 


= dynamic value of p on the p-y curve at depth x 
= corresponding reaction on the static p-y curve at 
depth x 

= frequency of loading, expressed in dimensionless 
terms a>r 0 /V s ; 

= circular frequency of loading equal to 2nf 
= lateral pile deflection relative to the soil at depth x 
= pile diameter (m) 
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a, (3, k and n are constants determined from curve fitting. 

In dynamic analysis, p-y curves are considered at 1 st , 4 th , 10 th 
and 30 th frequencies in soil-foundation interaction model 
with Chauk earthquake ground motion in longitudinal, 
transverse and vertical directions. 
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Fig. 6 Chauk Earthquake Motions 

IV. SOIL-FOUNDATION INTERACTION 

In the 3D soil modelling, the stiffness of the elastic 
surrounding soil is well accounted. The size of the soil model 
is 122m x 122m x 48m. The layered soil conditions are 
modelled using design parameters presented in Table IV and 
possion's ratio for sand is considered as 0.3 for all layers. 
The concrete pile tube is considered to behave linearly 
elastic and modelled as a cylindrical structure. The material 
properties of pile group are as follows: E=25GPa, y=23.58 
kN/m 3 and v=0.2. 

The interaction between the sand and the pile was modelled 
by defining tangential and normal contact behaviour in the 
finite element model. The tangential contact between the 
two surfaces was defined using a friction coefficient (tan 
(2(|)/3). The normal contact behavior of stiffness between the 
pile and soil was determined from p-y curves. The bottom of 
the pile and soil model was fixed and the exterior surface of 
the soil was constrained with multi-point constraints (MPC) 
constraints. 
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(d) 

Fig.8 Soil-foundation interaction (a) static, and dynamic 
analysis with Chauk earthquake in (b) Longitudinal, (c) 
Transverse and (d) Vertical direction with 0.12g of Chauk 
Earthquake 

In dynamic analysis in soil-foundation interaction, 
displacement in transverse direction of bridge is more than 
longitudinal and vertical direction of bridge in seismic 
excitations with Chauk earthquake ground motions. So, 
transverse direction of ground motion is considered with 
different frequencies. 


Fig. 7 Soil-foundation interaction finite element model 
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Fig.9 Soil-foundation interaction with earthquake 
excitation (a) 0.2g, and (b) 0.3g, (c) 0.4g and (d) 0.5g of 
Chauk Earthquake 

A. Analysis Results of Soil-Foundation Interaction in 
Static and Dynamic Conditions 

In static and dynamic conditions, maximum horizontal 
displacements are found in transverse direction of pier. 

Displacement in Horizontal and Vertical 


Pile 

(H) 

Pile 

(V) 

Cap 

(H) 

Cap 

(V) 

Soil 

(H) 



Different Motions 


Fig. 10 Comparison of static and dynamic analysis with 
different mode shapes of displacement of pile group and 
soil 
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Fig. 11 Comparison of static and dynamic analysis with 
different mode shapes (a) stresses and (b) strains 


vertical and horizontal displacements at the pile tip are 3mm 
and 2.94mm are found at 0.12g. After that, maximum ground 
acceleration of 0.5g is 14.5mm and 2.94mm in horizontal 
and vertical displacement of pile tip. Maximum shear stress 
and strain are found out the base of the pile cap. 

V. RESULTS AND DISCUSSIONS 

According to the soil investigation results presented Table II, 
the cohesion of soil, V for every layer are zero and the 
internal friction angle of soil are ranging from 19.5° to 30.5°. 
And then, soil resistance and deflection are determined for 
p-y curves. The p-y curves are developed for longitudinal 
and transverse directions of the foundation under studied 
pier. Due to the pile group configuration, the resistances of 
soil in transverse direction are less than longitudinal 
direction shown in Fig. 5. After developing the p-y curves, 
static and dynamic response of the soil-foundation 
interaction are carried out using ABAQUS. The analysis 
results are provided in Figs. 8 and 9. As shown in Fig. 9, the 
maximum displacements of the foundation of the pile cap in 
static and dynamic analysis are 0.2 mm and 14.4mm in 
horizontal direction of static and dynamic. And then, the 
maximum deflection of 0.14 and 14.5 mm are found out at 
the pile tip in longitudinal direction. On the other hand, the 
maximum vertical displacement of 2.28 mm and 3mm are 
found out at the pile tip. Due to the analysis results, the 
maximum principle stress and strain are found out at the 
base of the pile cap. The values of maximum principle stress 
and strain are 106658 N/m 2 and 4.8xl0 6 in static. And then, 
130900 N/m 2 and 5.84xl0 6 in dynamic of0.5g are found out 
respectively. Maximum and minimum principle stress and 
strain are found at the same region of the base of the pile cap 
and they are relating each other. 

VI. CONCLUSION 

This study focuses on the behaviour of soil-foundation 
interaction under a bridge pier, and the selected case study 
is pier number 5 (PR5) under the bridge of the middle span 
in Pakokku Bridge. According to the analysis results, the 
deflection of longitudinal direction is larger than transverse 
direction and the maximum deflection of 0.14 mm and 14.5 
mm are found out at the pile tip in longitudinal direction of 
static and dynamic condition of 0.5g. The maximum vertical 
displacement or settlement of 2.28mm and 3 mm are found 
out at the pile tip. From AASHTO Standard Specifications for 
Highway Bridges, 5'/L shall be limited to 0.005 for simple 
span bridges and 0.004 for continuous span bridges in 
vertical displacement and maximum horizontal 
displacement is 1 in [3]. The values of maximum principle 
stress and strain are 106658 N/m 2 and 4.8xl0 6 in static. 
And then, 130900 N/m 2 and 5.84xl0 6 in dynamic of 0.5g are 
found out respectively. In dynamic analysis, stress and strain 
in first mode are greater than the other mode shape. Finally 
it is found that the soil-foundation interaction under a bridge 
pier presented in this study is reliable and reasonable with 
the limitation of AASHTO Standard Specifications for 
Highway Bridges. 


According to analysis results with different mode shapes, the 
result of the first mode shape is approximately same as the 
other mode shapes. So, first mode shape is considered all of 
ground motions of Chauk earthquake. According to the 
analysis results in static condition, it is found that the 
vertical and horizontal displacements at the pile tip are 
2.28mm and 0.14mm respectively. In dynamic condition, the 
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